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We propose the model of a short-lived very powerful source of high 



^ ! Abstract 

^vq , energy neutrinos. It is formed as a result of the dynamical evolution 

CN | of a galactic nucleus prior to its collapse into a massive black hole 

and formation of high-luminosity AGN. This stage can be referred to 
as "pre- AGN". A dense central stellar cluster in the galactic nucleus 
on the late stage of evolution consists of compact stars (neutron stars 
Qh' and stellar mass black holes). This cluster is sunk deep into massive 

q . gas envelope produced by destructive collisions of a primary stellar 

f-j | population. Frequent collisions of neutron stars in a central stellar 

cluster are accompanied by the generation of ultrarelativistic fireballs 
and shock waves. These repeating fireballs result in a formation of the 
expanding rarefied cavity inside the envelope. The charged particles 
^ ' are effectively accelerated in the cavity and, due to pp-collisions in the 

gas envelope, they produce high energy neutrinos. All high energy 
particles, except neutrinos, are absorbed in the thick envelope. Du- 
ration of this pre-AGN phase is ~ 10 yr, the number of the sources 
can be ~ 10 per cosmological horizon. High energy neutrino signal 
can be detected by underground neutrino telescope with effective area 
S ~ 1 km 2 . 

PACS: 95.55.Vj, 95.85.Ry, 98.62.Js, 
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1 Introduction: Hidden Sources 



High energy (HE) neutrino radiation from astrophysical sources is accompa- 
nied by other types of radiation, most notably by the HE gamma-radiation. 
These HE gamma-radiation can be used to put upper limit on the neutrino 
flux emitted from a source. For example, if neutrinos are produced due to in- 
teraction of HE protons with low energy photons in extragalactic space or in 
the sources transparent for gamma-radiation, the upper limit on diffuse neu- 
trino flux I V {E) can be derived from e-m cascade radiation. This radiation is 
produced due to collisions with photons of microwave radiation j bb , such as 
7 + 7^ — > e + + e~, e + %b "~ > e ' + 7' etc. These cascade processes transfer 
the energy density released in high energy photons co> 7 into energy density 
of the remnant cascade photons u cas . These photons get into the observed 
energy range 100 MeV— 10 GeV and their energy density is limited by recent 
EGRET observations U as uo cas < 2 ■ 10~ 6 eV/cm 3 . Introducing the energy 
density for neutrinos with individual energies higher than E, u u (> E), it is 
easy to obtain the following chain of inequalities (reading from left to write) 

00 00 

u cas > u u {> E) = — f EI v {E)dE >—E( I v {E)dE = —EI U {> E). (1) 
c J c J c 

E E 

Now the upper limit on the integral HE neutrino flux can be written down 

as 

/„(> E) < — ^ = 4.8 ■ lO 3 ^ 1 cm-VV- 1 . (2) 
c E 

However, there can be sources, where accompanying electromagnetic radi- 
ation, such as gamma and X-rays, are absorbed. They are called "hidden 
sources" 0. Several models of hidden sources were discussed in the litera- 
ture. 

• Young SN Shell 0] during time t„ ~ 10 3 — 10 4 s are opaque for all 
radiation, but neutrinos. 

• Thome- Zytkow Star the binary with a pulsar submerged into a 
red giant, can emit HE neutrinos while all kinds of e-m radiation are 
absorbed by the red giant component. 

• Cocooned Massive Black Hole (MBH) in AGN @ is an example of 
AGN as hidden source: e-m radiation is absorbed in a cocoon around 
the massive black hole. 
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• AGN with Standing Shock in the vicinity of a MBH || can produce 
large flux of HE neutrinos with relatively weak X-ray radiation. 

In this paper we propose a new model of the hidden source which can op- 
erate in a galactic nucleus at pre-AGN phase, i.e. prior to MBH formation in 
it. The MBH in AGN is formed through the dynamical evolution of a central 
stellar cluster resulting in a secular contraction of the cluster and its final 
collapse f7|, ||. The first stage of this evolution is accompanied by collisions 
and destruction of normal stars in the evolving cluster, when virial velocities 
of constituent stars become high enough. The compact stars (neutron stars 
and black holes) survive this stage and their population continue to con- 
tract, being surrounded by the massive envelope composed of the gas from 
destroyed normal stars. Pre-AGN phase corresponds to a near collapsing 
central cluster of compact stars in the galactic nucleus. Repeating fireballs 
after continuous collisions of compact stars in this very dense cluster result in 
the formation of a rarefied cavity in the massive gas envelope. Particles ac- 
celerated in this cavity interact with the gas in the envelope and produce HE 
neutrinos. Accompanying gamma-radiation can be fully absorbed in the case 
of thick envelope (matter depth X env ~ 10 4 g/cm 2 ). The proposed source is 
short-lived (lifetime t s ~ 10 years) and very powerful: neutrino luminosity 
exceeds the Eddington limit for e-m radiation. 



2 The Model 

We consider in the following the basic features of the formation of a short- 
lived extremely powerful hidden source of HE neutrinos in the process of 
dynamical evolution of the central stellar cluster in a typical galactic nucleus. 



2.1 Dynamical Evolution of Galactic Nucleus 

The dynamical evolution of dense central stellar clusters in the galactic nuclei 
is accompanied by a secular growth of the velocity dispersion of constituent 
stars v or, equivalently, by the growth of the central gravitational potential. 
This process is terminated by the formation of the MBH when the velocity 
dispersion of stars grows up to the light speed (see for a review e.g. || and 
references therein). On its way to a MBH formation the dense galactic nuclei 
inevitably proceed through the stellar collision phase of evolution p|-|[13|, 
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when most normal stars in the cluster are disrupted in mutual collisions. The 
necessary condition for the collisional destruction of normal stars with mass 
m* and radius r* in the cluster of identical stars with a velocity dispersion v 
is 

v > v p , (3) 

where 

f2Gm,\ 1 ' 2 2 ( m.V' 2 f >: Y 1 ' 2 ,i 

is an escape (parabolic) velocity from the surface of a constituent normal star. 
The kinetic energy of colliding star is greater in general than its gravitational 
bound energy under the inequality (|3|). If v > v p , the normal stars are 
eventually disrupted in mutual collisions or in collisions with the extremely 
compact stellar remnants, i.e. with neutron stars (NSs) or stellar mass black 
holes. Only these compact stellar remnants will survive through the stellar- 
destruction phase of evolution (v = v p ) and form the self-gravitating core. 
We shall refer for simplicity to this core as to the NS cluster. Meanwhile the 
remnants of disrupted normal stars form a gravitationally bound massive 
gas envelope in which the NS cluster is submerged. The virial radius of this 
envelope is 

GM env lM env ( Tn* \ ~^ / \ 

R ™ = ^H = 4^7 r * * °- 56Ms \W & ) \r & ) pc ' (5) 

where M env = 10 8 M 8 M Q is a corresponding mass of the envelope. The gas 
from disrupted normal stars composes the major part of the progenitor cen- 
tral stellar cluster in the galactic nucleus. So the natural range for the total 
mass of the envelope is the same as the typical range for the mass of a central 
stellar cluster in the galactic nucleus, M env = 10 7 — 1O 8 M . The envelope 
radius R env is given by the virial radius of a central cluster in the galactic nu- 
cleus at the moment of evolution corresponding to normal stars destructions, 
i. e, v — v p . The mean number density of gas in the envelope is 

n env = A_L^ ~ 5.4 . i 9M 8 " 2 (—X(-) 3 cm- 3 , (6) 



where m p is a proton mass. A column density of the envelope is 

M 



■p 

X env = m p n env R env ~ 1.6 ■ 10 4 M 8 _1 I —A I ( ^- I gem -2 . (7) 
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Such an envelope completely absorbs electromagnetic radiation and HE par- 
ticles outgoing from the interior, except neutrinos and gravitational waves. 
A column density becomes less for more massive envelopes. 

We assume that energy release due to star collisions supports the gas 
in the cluster in (quasi-) dynamical equilibrium. It implies the equilibrium 
temperature T eg of the gas, 

^GMem,^ , K (8) 

The thermal velocity of gas particles at the equilibrium temperature is of 

order of an escape velocity from the surface of a normal star. 

If T ^> T eg the gas outflows from the cluster with the sound speed 

/ kT\ 1/2 (T\ 1/2 
Vs = 27 — « 600 — km/s, (9) 



nipj \T t 



where 7 is adiabatic index (7 = 5/3 for hydrogen). 
If T <C T eq gas collapses to the core. 



2.2 Dense Cluster of Stellar Remnants 

As was discussed above, the dense NS cluster survives inside the massive 
envelope of the post-stellar-destruction galactic nucleus. The total mass of 
this cluster is ~ 1 — 10 % of the total mass of a progenitor galactic nucleus 



P|-|12|| and so of the massive envelope, i. e. M ~ 0.01 — 0.1M env . We will use 
the term 'evolved galactic nucleus' for this cluster of NSs assuming that (i) 
v > v p and (ii) the (two-body) relaxation time in the cluster is much less than 
the age of the host galaxy. Under the last condition the cluster has enough 
time for an essential dynamical evolution. For example the relaxation time of 
stars inside a central parsec of the Milky Way galaxy is t r ~ 10 7 — 10 8 years. 
A further dynamical evolution of the evolved cluster is terminated by the 
dynamical collapse to a MBH. 

We consider in the following an evolved central cluster of NSs with iden- 
tical masses m = 1.4M . This evolved cluster of NSs is sunk deep into the 
massive gas envelope remaining after the previous evolution epoch of a typ- 
ical normal galactic nucleus. Let N = M/m = lO 6 ^ is a total number of 
NSs stars in the cluster. The virial radius of this cluster is: 

R=^^ = - -) Nr q ~ 1.0- 10 13 iV 6 (w/0.1c)- 2 cm, (10) 
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where r g = 2Gm/c 2 is a gravitational radius of NS. For N ~ 10 6 and v ~ 0.1c 
one has nearly collapsing cluster with the virial size of ~ 1 AU. The charac- 
teristic times are (i) the dynamical time td yn — R/v — (l/4:)N(c/v) 3 r g /c ~ 
0.95Nq(v /0.1c) -3 hour and (ii) the evolution (two-body relaxation) time of 
the NS cluster t rd ~ 0.1(N/ In N)t dyn ~ 19iVf (f/O.lc)- 3 years. In general 
trei 3> tdyn, if ./V 3> 1. This evolution time determines the duration of an 
active phase for the considered below hidden source, as t s ~ t re[ ~ 10 years. 



2.3 Fireballs in Cluster 

The most important feature of our model is a secular growing rate of acciden- 
tal NS collisions in the evolving cluster, accompanied by large energy release. 
The corresponding rate of NS collisions in the cluster (with the gravitational 
radiation losses taken into account) is JTJJ-[|njj 



v\ 17 / 7 c rf v \ 31/7 1 c 



) R =6K)yi \-c) Nr 



fi 



~ 4.4- itfN^iv/U.lcf 1 ' 7 yr -1 . (11) 

The time between two successive NS collisions is 

. . 1 /c\ 10 / 7 1 /c\ 31 /7 r 

~ 7.3 ■ 10 3 A^ 6 (i;/0.1c)" 31/7 s. (12) 

Note that a number of NS collisions, given by N c t s , comprises only a small 
fraction, about 1%, of a total number of NSs in the cluster to the time of the 
onset of dynamical collapse of the whole cluster into a MBH. 

Merging of two NSs in collision is similar to the process of a tight binary 
merging: the NSs are approaching to each other by spiralling down due to 
the gravitational waves radiation and then coalesce by producing an ultrarel- 



ativistic photon-lepton fireball [0-0], which we assume to be spherically 
symmetric. 

The energy of one fireball is E = E 52 10 52 ergs, and the total energy 
release in the form of fireballs during lifetime of the hidden source t s ~ 10 yr 
is 

E tot ~ N c E t s ~ 4 ■ 10 56 ergs, (13) 
where N r is the NS collision rate. 
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The physics of fireballs is extensively elaborated especially in recent years 
for the modeling of cosmological gamma-ray bursts (GRBs) (for review see 
e. g. |22| and references therein). The newborn fireball expands with rela- 
tivistic velocity, corresponding to the Lorentz factor Tf ^> 1. The relevant 
parameter of a fireball is the total baryonic mass 

M = E /r]c 2 ~ 5.6 ■ KT 6 £ 52 % ^ (14) 

where baryon- loading mass parameter rj = lO 3 ^. The maximum possible 
Lorentz factor of expanding fireball is Tf = rj+1 during the matter-dominated 



phase of fireball expansion [H], |2(J. During the initial phase of expansion, 
starting from the radius of the 'inner engine' Ro ~ 10 6 — 10 7 cm, the fireball 
Lorentz factor increases as T oc r, until it is saturated at the maximum 
value Tf = rj ^ 1 at the radius R v = RqI] (see e. g. P2fl). Internal shocks 
will take place around R s h = Rot] 2 , if the fireball is inhomogeneous and the 
velocity is not a monotonic function of radius, e.g. due to the considerable 



emission fluctuations of the inner engine p3|-||25||. The fireball expands with 
the constant Lorentz factor r = r\ at R > R v until it sweeps up the mass 
Mq/tj of ambient gas and looses half of its initial momentum. At this moment 



(R = _R 7 ) the deceleration stage starts pl|-p7|. 

Interaction of the fireball with an ambient gas determines the length of its 
relativistic expansion. In our case the fireball propagates through the massive 
envelope with a mean gas number density n env = p e nv/m P = n$10 9 cm" 3 as 
it follows from Eq. (|). Fireball expands with r > 1 up to the distance 
determined by the Sedov length 

l S = f f A?) - 12 • 10 1 V /3 i4 /3 cm. (15) 

y47T p env C ) 

Fireball becomes mildly relativistic at radius r = 1$ due to sweeping up the 
gas from the envelope with the mass MqT]. The radius r = Is is the end 
point of the ultrarelativistic fireball expansion phase. Far beyond the Sedov 
length radius (r ^> Is) there is the non-relativistic Newtonian shock driven by 
the decelerated fireball. Its radius R(t) obeys the Sedov-Taylor self-similar 
solution f28|| , with R(t) = (E t 2 / Penv) 1 ^ ■ The corresponding shock expansion 



velocity is u = (2/5) [l s /R{t)f /2 c < c. 
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2.4 Cavity and Shocks 



We show here that relativistic fireballs from a dense central cluster of NSs 
produce the dynamically supported rarefied cavity deep inside the massive 
gaseous envelope. 

The first fireball sweeps out the gas from the envelope producing the 
cavity with a radius lg. This cavity expands due to the next fireballs, which 
propagate first in a rarefied cavity and then hit the boundary pushing it 
further. 

Each fireball hitting the dense envelope is preceded by a shock. Propa- 
gating through the envelope, the shock sweeps up the gas ahead of it and 
gradually decelerates. The swept out gas forms a thin shell with a density 
profile given by Sedov self-similar solution. The next fireball hits this thin 
shell when it is decelerated down to the non-relativistic velocity. Moving in 
the envelope, the shell accumulates more gas, retaining the same density pro- 
file, and then it is hit by the next fireball again. After a number of collisions 
the shell becomes massive, and the successive hitting fireballs do not change 
its velocity appreciably. In this regime one can consider the propagation of 
massive non-relativistic thin shell with a shock (density discontinuity) ahead 
of it. The shock speed v s h is connected with a velocity v g of gas behind it as 
v sh = (7 + l)fg/2, where 7 is an adiabatic index. The density perturbation 
in the envelope propagates as a shock until v s h remains higher than sound 
speed v s . In the considered case, the shock dissipates in the middle of the 
envelope. For the Sedov solution the shock velocity changes with distance r 
as 



where as is the constant of self-similar Sedov solution; when radiative pres- 
sure dominates as = 0.894. The other quantities in Eq. ( |i~6|) are E sh = 
(1/2) Et t = 2 ■ 10 56 erg/s is the energy of the shock, which includes kinetic 
and thermal energy (the half of a total energy is transformed to particles 
accelerated in the cavity), and p env is a density of the envelope given by 
Eq. (|). From Eqs. fll6|) and flj^) it follows that v s h > v s holds at distance 
r < 1 ■ 10 18 cm ~ 0.6R env , i. e. that shock does not reach the outer surface 
of the envelope. In fact the latter conclusion follows already from energy 
conservation. The gravitational energy of the envelope Vq is 




(16) 




env 



> 9.2 • 10 &b ergs, 



(17) 
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where k depends on a radial profile of the gas density in the envelope p(r), and 
changes from 3/5 to 1. This energy is higher than a total injected energy 
Etot ~ 4 ■ 10 56 ergs, and thus the system remains gravitationally bound. 
When a shock reaches the boundary of the envelope, the gas distribution 
changes there. It has the form of a thin shell with gravitational energy 
V e = GM^ nv /2R env . To provide the exit of the shock to the surface of the 
envelope, a total energy release must satisfy the relation E tot > V — V e , 
where the minimum value of Vo — V e is 1.5 • 10 56 ergs for k = 3/5. Actually 
Etot must be higher because (i) part of the injected energy goes to heat, 
(ii) more realistic value of k = 1, and (iii) the shell has a non-zero velocity, 
when the shock disappears, with a gravitational braking taken into account. 
We conclude thus that for E tot ~ 4 • 10 56 ergs, shock dissipates inside the 
envelope. 

The cavity radius grows with time. For the stage, when a shell moves 
non-relativistically, the cavity radius calculated as a distance to the shell in 
the Sedov solution is 



jw*)= (-^y v. (is) 



At the end of a phase of the hidden source activity, t s ~ 10 yr, the radius of 
the cavity reaches ~ 3 • 10 17 cm, remaining thus much less than R env . 

The cavity is filled by direct and reverse relativistic shocks from fireballs. 
Reverse shocks are produced by decelerated fireballs, most notably when they 
hit the boundary of the cavity. The expanding fireballs inside the cavity have 



a shape of thin shells 29 and are separated by distance 



R c = ct e ~ 2.2 ■ 10 14 A" 6 (w/0.1c)- 31/7 cm. (19) 

The gas between two fireballs is swept up by the preceding one. A total 
number of fireballs existing in the cavity simultaneously is Nf ~ R cav /R c ^> 
1, and this number grows with time as t 2 ^ . 

Shocks generated by repeating fireballs are ultrarelativistic inside the cav- 
ity and mildly relativistic in the envelope near the inner boundary. Colli- 
sions of multiple shocks in the cavity, as well as inside the fireballs, produce 
strongly turbulized medium favorable for generation of magnetic fields and 
particle acceleration. 
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3 High Energy Particles in Cavity 



There are three regions where acceleration of particles take place. 

(i) NS cluster, where fireballs collide, producing the turbulent medium 
with large magnetic field. This region has a small size of order of virial radius 
of the cluster R ~ 10 13 cm, and we neglect its contribution to production of 
accelerated particles. 

(ii) The region at the boundary between the cavity and the envelope. 
During the active period of a hidden source, t s ~ 10 yr, the fireballs hit this 
region, heating and turbulizing it. The large magnetic equipartition field is 
created here. This boundary region has density p ~ p enV) the radius R ~ R cav 
and the width A < 0.1-R ca „. 

(iii) Most of the cavity volume is occupied by fireballs, separated by 
distance R c . Due to collisions of internal shocks, the gas in a fireball is 



turbulized and equipartition magnetic field is generated [22], |30 . 

In all three cases the Fermi II acceleration mechanism operates. For all 
three sites we assume existence of equipartition magnetic field, induced by 
turbulence and dynamo mechanism: 

8^ ~ — ' (20) 

where p and u t are the gas density and velocity of turbulent motions in the 
gas, respectively. Since turbulence is caused by shocks, the shock spectrum 
of turbulence Ff. ~ k~ 2 is valid, where A; is a wave number. Assuming 
equipartition magnetic field on each scale I ~ Hf ~ kF^, one obtains 

the distribution of magnetic fields over the scales as 

Ht/Ho = (Z/Z ) 1/2 , (21) 

where Iq is a maximum scale with the coherent field Hq there. 
The maximum energy of accelerated particles is given by 

E m ax ~ eH l (22) 

with an acceleration time 

< 23 > 

For the turbulent shell at the boundary between the cavity and envelope, 
assuming mildly relativistic turbulence u t ~ c and p ~ p env , we obtain 
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H eq = 4 • 10 3 G. The maximum acceleration energy is E max = 2 • 10 21 eV, if 
the coherent length of magnetic field Iq is given by the Sedov length 1$, and 
the acceleration time is t acc = 4 ■ lO^i^Tig 1//3 s. The typical time of energy 
losses, determined by pp-collisions, is much longer than t acc , and does not 
prevent acceleration to E max given above: 

t fldE\~ 1 _ 1 _ 2 10 6 n -i s (24) 

where f p ~ 0.5 is the fraction of energy lost by HE proton in one collision, 
(jpp is a cross-section of pp-interaction, and n env is the gas number density in 
the boundary turbulent shell. 

The turbulence in the fireball is produced by collisions of internal shells, 
and a natural scale for coherence length 1$ is a width of the internal shell in the 
local frame 5'. The maximum energy in the laboratory frame E max ~ eH' eq 5, 
where 5 is the corresponding width in the laboratory frame. Since H' oc 1/R 
and 5 oc R, the maximum energy does not change with time, and can be 
estimated as in Ref. |3(J, E 

max ~ 3 • 10 2 ^ eV. Note that in our case a fireball 



propagates in the very low-density medium. 

A gas left in the cavity by preceding fireball, as well as high energy 
particles escaping from it, are accelerated by the next fireball by a factor 
r 2 at each collision. This r 2 -mechanism of acceleration works only in pre- 
hydrodynamic regime of fireball expansion, after reaching the hydrodynamic 



stage, T -mechanism ceases [31 



We conclude, thus, that both efficiency and maximum acceleration en- 
ergy are very high. We assume that a fireball transfers half of its energy to 
accelerated particles. 



4 Neutrino Production and Detection 

Particles accelerated in the cavity interact with the gas in the envelope pro- 
ducing high energy neutrino flux. We assume that about half of the total 
power of the source L tot is converted into energy of accelerated particles 
L p ~ 7 ■ 10 47 erg/s. As estimated in Section 2.1, the column density of 
the envelope varies from X env ~ 10 2 g/cm 2 (for very heavy envelope) up to 
X env ~ 10 4 g/cm 2 (for the envelope with mass M ~ 1O 8 M . Taking into 
account the magnetic field, one concludes that accelerated protons loose in 
the envelope a large fraction of their energy. The charged pions, produced in 
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pp-collisions, with Lorentz factors up to T c ~ 1/ ((^ W N n envC T w) ~ 4 • lO 13 ^ 1 
freely decay in the envelope (here a n N ~ 3 • 10~ 26 cm 2 is 7riV-cross-section, r n 
is the lifetime of charged pion, and n env = I0 9 n$ cm -3 is the number density 
of gas in the envelope). We assume E~ 2 spectrum of accelerated protons 

Q P (E) = (25) 

where ( = \a{E max / E min ) ~ 20—30. About half of its energy protons transfer 
to high energy neutrinos through decays of pions, L v ~ (2/3)(3/4)L p , and 
thus the production rate of + neutrinos is 

Q^(>E) = -^. (26) 

Crossing the Earth, these neutrinos create deep underground the equilibrium 
flux of muons, which can be calculated as |32| : 

where the normalization cross-section oo = 1 • 10~ 34 cm 2 , A^^ = 6 • 10 23 is 
the Avogadro number, 6 M = 4 • 10 -6 cm 2 /g is the rate of muon energy losses, 
Y^(E) is the integral muon moment of u^N interaction (see e. g. [0, |32]| ). 
The most effective energy of muon detection is > 1 TeV . The rate of 
muon events in the underground detector with effective area S at distance r 
from the source is given by 

Thus, we expect about 10 muons per year from the source at distance 
10 3 Mpc. 



5 Accompanying Radiation 

We shall consider below HE gamma-ray radiation produced by accelerated 
particles and thermalized infrared radiation from the envelope. As far as HE 
gamma-ray radiation is concerned, there will be considered two cases: (i) 
thin envelope with X env ~ 10 2 g/cm 2 and (ii) thick envelope with X env ~ 
10 4 g/cm 2 . In the latter case HE gamma-ray radiation is absorbed. 
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5.1 Gamma- Ray Radiation 



Apart from high energy neutrinos, the discussed source can emit HE gamma- 
radiation through 7T° — > 27 decays and synchrotron radiation of the electrons. 
In case of the thick envelope with X env ~ 10 4 g/ cm 2 most of HE photons are 
absorbed in the envelope (characteristic length of absorption is the radiation 
length X rad « 60 g/cm 2 ). In the case of the thin envelope, X env ~ 100 g/cm 2 , 
HE gamma-radiation emerges from the source. Production rate of the syn- 
chrotron photons can be readily calculated as 

dE 

dQ syn = -^Q e {> E e ), (29) 

where E e and E 1 are the energies of electron and of emitted photon, respec- 
tively. Using Q e (> E e ) = L e /(r)E e ) and E 7 = k syn (H)E^, where k syn is the 
coefficient of the synchrotron production, one obtains 

Qsyn(E 7 ) = y^^- (30) 



Note, that the production rate given by Eq. ( p0|) does not depend on magnetic 
field. Adding the contribution from n° — > 27 decays, one obtains 

L, 



= 12^' (31) 
and the flux at E~ > 1 GeV at the distance to the source r = 1 • 10 3 Mpc is 



F 7 (> E 7 ) 




-2 

-2 -1 



10 47 erg s- 1 J \10 3 Mpc / 
i. e. the source is detectable by EGRET. 



cm' 2 s" 1 , (32) 



5.2 Infrared and Optical Radiation 

Hitting the envelope, fireballs dissipate part of its kinetic energy in the enve- 
lope in the form of low-energy e-m radiation. This radiation is thermalized 
in the optically thick envelope and then re-emitted in the form of black- 
body radiation from the surface of the envelope. It appears much more later 
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than HE neutrino and gamma-radiation. The thermalized radiation diffuses 
through the envelope with a diffusion coefficient D ~ cZ^/, where a diffusion 
length is Idif = l/(crTn env ) and ot is the Thompson cross-section. A mean 
time of the radiation diffusion through the envelope of radius R env is 

t d ~ ^ ~ 1 • 10 4 yr, (33) 

independently of the envelope mass. This diffusion time is to be compared 
with a duration of an active phase t s ~ 10 years and with a time of flight 
Renv/c ~ 2 years. 

Since the produced burst is very short, t s ~ 10 yr, the arrival times of 
photons to the surface of the envelope have a distribution with a dispersion 
a ~ td- An average surface black body luminosity is then 

L bb ~ E tot /t d ~ 1 • 10 45 erg/s, (34) 

with a peak luminosity being somewhat higher. The temperature of this 
radiation corresponds to IR range 

(J \ 1/4 
T bb = I— »_ ^8.4-10 2 K. (35) 

\^Rl nv (JSBj 

Thus, in ~ 10 4 years after neutrino burst, the hidden source will be seen in 
the sky as a luminous IR source. Since we consider a model of production of 
high-luminosity AGN, the object is typically expected to be at high redshift 
z. Its visible magnitude is 

where H = lOO/i km/s Mpc is the Hubble constant and the flux f = 
2.48 • 10" 5 erg/cm 2 s. For L bb = 1 • 10 45 erg/s, z = 3 and h = 0.6, the visible 
magnitude of the IR source is m = 22.7. Such a faint source is not easy 
identify in e.g. IRAS catalogue as a powerful source, because for redshift 
determination it is necessary to detect the optical lines from the host galaxy, 
which are very weak at assumed redshift z — 3. The non-thermal optical 
radiation can be also produced due to HE proton-induced pion decays in the 
outer part of the envelope, but its luminosity is very small. Most probably 
such source will be classified as one of the numerous non- identified weak IR 
source. 
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5.3 Duration of activity and the number of sources 

As was indicated in Section 2.2 the duration of the active phase t s is deter- 
mined by relaxation time of the NS cluster: t s ~ t re i ~ 10 — 20 yr. This 
stage appears only once during the lifetime of a galaxy, prior to the MBH 
formation. If to assume that a galactic nucleus turns after it into AGN, the 
total number of hidden sources in the Universe can be estimated as 

4 

N H s ~ -K(3ct ) 3 n AGN t s /t AGN , (37) 

where |"7r(3cto) 3 is the cosmological volume inside the horizon cto, ^agn is the 
number density of AGNs and t A GN is the AGN lifetime. The value t s /tAGN 
gives a probability for AGN to be observed at the stage of the hidden source, 
if to include this short stage (t s ~ 10 yr) in the much longer (t^Gv) AGN 
stage and consider (for the aim of estimate) the hidden source stage as the 
accidental one in the AGN history. The estimates for uagn and tAGN taken 
for different populations of AGNs result in N HS ~ 10 — 100. 



6 Conclusions 

Dynamical evolution of the central stellar cluster in the galactic nucleus re- 
sults in the stellar destruction of the constituent normal stars and in the 
production of massive gas envelope. The surviving subsystem of NSs sub- 
merges deep into this envelope. The fast repeating fireballs caused by NS 
collisions in the central stellar cluster produce the rarefied cavity inside the 
massive envelope. Colliding shocks generate the turbulence inside the fire- 
balls and in the cavity, and particles are accelerated by Fermi II mechanism. 
These particles are then re-accelerated by r 2 -mechanism in collisions with 
relativistic shocks and fireballs. 

All high energy particles, except neutrinos, can be completely absorbed 
in the thick envelope. In this case the considered source is an example of a 
powerful hidden source of HE neutrinos. 

Prediction of high energy gamma-ray flux depends on the thickness of 
envelope. In case of the thick envelope, X ra d ~ 10 4 g/cm 2 , HE gamma- 
radiation is absorbed. When an envelope is thin, X rad ~ 10 2 g/cm 2 , gamma- 
ray radiation from 7r° — > 27 decays and from synchrotron radiation of the 
secondary electrons can be observed by EGRET and marginally by Whipple 
detector at -EL > 1 TeV. 
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In all cases the thickness of the envelope is much larger than the Thomp- 
son thickness (xt ~ 3 g/cm 2 ), and this condition provides the absorption 
and X-rays and low energy gamma-rays. 

A hidden source is to be seen as a bright IR source but, due to slow diffu- 
sion through envelope, this radiation appears ~ 10 4 years after the phase of 
neutrino activity During the period of neutrino activity the IR luminosity 
is the same as before it. A considered source is a precursor of most powerful 
AGN, and therefore most of these sources are expected to be at the same 
redshifts as AGNs. The luminosity Lm ~ 10 45 — 10 46 erg/s is not unusual 
for powerful IR sources from IRAS catalogue. The maximum observed lumi- 
nosity exceeds 1 • 10 48 erg/s [33|, and there are many sources with luminosity 



10 45 — 10 46 erg/s [Q]. Moreover, for most of the hidden sources the distance 
cannot be determined, and thus they fall into category of faint non-identified 
IR sources. 

Later these hidden sources turn into usual powerful AGNs, and thus the 
number of hidden sources is restricted by the total number of these AGNs. 

In our model the shock is fully absorbed in the envelope. Since the total 
energy release E tot is less than gravitational energy of the envelope E grav ~ 
GMl nv j R env , the system remains gravitationally bound, and in the end the 
envelope will collapse into black hole or accretion disc. 

The expected duration of neutrino activity for a hidden source is ~ 10 yr, 
and the total number of hidden sources in the horizon volume ranges from a 
few up to ~ 100, within uncertainties of the estimates. 

Underground neutrino detector with an effective area S ~ 1 km 2 will 
observe ~ 10 muons per year with energies > 1 TeV from this hidden 
source. 
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